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CONVERGENCE TO EQUILIBRIUM FOR THE KINETIC 
FOKKER-PLANCK EQUATION ON THE TORUS 


HELGE DIETERT, JOSEPHINE EVANS, AND THOMAS HOLDING 

Abstract. We study convergence to equilibrium for the kinetic Fokker-Planck 
equation on the torus. Solving the stochastic differential equation, we show 
exponential convergence in the Monge-Kantorovich-Wasserstein W 2 distance. 
Finally, we investigate if such a coupling can be obtained by a co-adapted 
coupling, and show that then the bound must depend on the square root of 
the initial distance. 


1. Introduction 


The kinetic Fokker-Planck equation, also known as the Kramers equation, is a 
basic model for the spreading of a solute due to interaction with the fluid back¬ 
ground. It is derived from Langevin dynamics, where the time scale of observation is 
much larger than the correlation time of the solute-fluid interactions (see e.g. IT511. 
In the context of fixed random scatters the similar linear Landau equation can rig¬ 
orously be derived in the weak coupling limit, see ||B] and references within. We 
focus on the case that the space variable is in the torus T = R/(27rLZ) of length 
2irL. The kinetic Fokker-Planck equation describes the law of a particle moving in 
the phase space Txl whose location in the phase space is ( X t , Vj) and evolves as 



( 1 ) 


where dWt is a standard white noise. The corresponding measure (i t on I x 1 
evolves as 


( 2 ) 


dtUt + vd x nt = d v [Xvut + yd v fj, t ], 


where this equation is considered in the weak sense. 

We expect that the measure /x* spreads out over time and eventually reaches the 
uniform measure which is the unique stationary state. The problem of convergence 
to equilibrium has been studied in different metrics before (see e.g. [51[TT]1 and 
forms a key example of hypocoercivity. 

To the best of the authors’ knowledge convergence in the Monge-Kantorovich- 
Wasserstein (MKW) distance W 2 has not been solved and is the object of this 
paper. 

The MKW distance comes from optimal transport and is defined as 
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where IIis the set of all couplings between fi and v. This metric is very useful 
as it allows to understand the Fokker-Planck equation as gradient flow [3], see m 
for a general review. 

In the spatially homogeneous case, i.e. only considering V t . this is an Orstein- 
Uhlenbeck process for which exponential convergence to equilibrium has been proven 
in the Wasserstein distance [I]. In a stochastic framework the convergence can be 
proved by coupling the noise and using the fact that the dependence on the ini¬ 
tial data decays over time, which in an analytic setting translates to a functional 
inequality for the time derivative showing that the evolution is a contraction semi¬ 
group. 

In the case where there is also a spatial variable, the same coupling approach 
works if the spatial variable evolves in a confining potential. However, in our case 
on the torus the spatial distance will not decay if we just couple the velocities. 

Solving the stochastic evolution, we are still able to show exponential decay of 
the distance between two solutions. 

Theorem 1. If fit and v t are two solutions to the kinetic Fokker-Planck equation 
ED, then we have 

< (e~ xt +ce- t/4A2L2 ) W 2 (/x 0 , v 0 ) 
for a constant c only depending on L. 

The key idea is that, after fixing the net effect of the velocity noise, the spatial 
variable has enough randomness left to allow such a coupling. This approach is 
not based on a functional inequality which is integrated over time and in fact the 
evolution is not a contraction semigroup. 

Theorem 2. There exists no 7 > 0 such that for all solutions fit and v t to the 
kinetic Fokker-Planck equation ED we have 

W 2 (fit,vt) < e -7t W 2 (/Ai, v o) 

for all t > 0. 

This shows that the generator is not coercive but only hypocoercive in W 2 . 

In probability theory a classical approach to such convergence results is the 
construction of a coupling [8]. For this, random variables (Xfi Vf ) are constructed 
for t £ R+ and * = 1,2 such that (X^V/) has law fi t and (X 2 ,V 2 ) has law v t . 
Then for t £ R + the coupling ((Xf,V t 1 ),(X 2 ,V 2 )) gives an upper bound of the 
MKW distance W 2 (/Ut, v t ). 

The stochastic differential equation ED motivates to look at couplings where 
(XI, Vf) are continuous Markov processes with initial distribution fiQ and vq, re¬ 
spectively, and whose transition semigroup is determined by ED- For such couplings 
we can consider a more restrictive class of couplings. 

Definition 3 (co-adapted coupling). The coupling ((Xf, Vf ), (Xf, V 2 )) is co-adapted 
if, fori = 1, 2, under the filtration T generated by the coupling ((Xf, V) 1 ), (X 2 , Vfi)), 
the process (XfiVf ) is a continuous Markov process whose transition semigroup is 
determined by ED- 

This is an important subclass of couplings, which contains many natural cou¬ 
plings, and an even more restrictive subclass is the class of Markovian couplings, 
where additionally the coupling itself is imposed to be Markovian. The existence 
and obtainable convergence behaviour under this restriction has already been stud¬ 
ied in different cases, e.g. HUME Note that the co-adapted coupling is equivalent 


CONVERGENCE FOR THE KINETIC FOKKER-PLANCK EQUATION 


3 


to the condition that the filtration generated by ( X \, Vf) is immersed in the filtra¬ 
tion generated by the coupling, which motivates Kendall 0 to call such couplings 
immersed couplings. 

By adapting the reflection/synchronisation coupling, we can still obtain expo¬ 
nential convergence but with a loss in dependence on the initial data. 


Theorem 4. Given initial distributions po and v$, there exists a co-adapted cou¬ 
pling {{X}, I/, 1 ), (X 2 , V 2 )) such that 

< (E [IX, 1 - X 2 | 2 + (K/ - V t 2 ) 2 ] ) 1/2 

< Cj3(t)(y/W2{p o, vq) + W 2 {po, vq)), 


where 


P(t) = 


-min(2A,l/(2A 2 L 2 ))t 

e-2A*(l +t) 


and C is a constant that depends only on A and L. 


4L 2 A 3 ^ 1 
4L 2 A 3 = 1 


Here we used the notation jX 1 — Xf\j to emphasis that this is the distance on 
the torus T. In fact the filtration generated by (X 1 , V 1 ) and (X 2 , V 2 ) agree which 
Kendall j5] calls an equi-filtration coupling. 


Remark 5. This achieves the same exponential decay rate as the non-Markovian 
argument, except for the case 4L 2 A 3 = 1, when the spatial and velocity decay rates 
coincide and we have an addition polynomial factor. 


In general the loss in the dependence is necessary. 

Theorem 6. Suppose there exists a function a : t— > M. + and a constant 7 > 0 

such that for all initial distributions po and v$ there exists a co-adapted coupling 
{{X},V?),{X 2 ,V 2 )) such that 

(E [IX, 1 - X 2 | 2 + (K, 1 - V t 2 ) 2 ]) 1/2 < a(W 2 (p 0 ,v 0 ))e-^. 

Then there exists a constant C such that for z G (0,7rL] we have the following lower 
bound on the dependence on the initial distance 

a(z) > Cyfz. 

The idea is to focus on a drift-corrected position on the torus, which evolves as 
a Brownian motion. By stopping the Brownian motion at a large distance we can 
then prove the claimed lower bound. 

This shows that a simple hypocoercivity argument on a Markovian coupling can¬ 
not work. Precisely, there cannot exist a semigroup P on the probability measures 
over (Txl) x2 , whose marginals behave like the solution of <[T|) and which satisfies 
H(P t {n)) < cH{n)e~^ for H 2 (n) = /[(X 1 -X 2 ) 2 + (K 1 -I/ 2 ) 2 ]d 7 r(X 1 , V 1 , X 2 , V 2 ). 
Otherwise, the Markov process associated to P would be a coupling contradicting 
Theorem 


1.1. Acknowledgements. The authors would like to thank Clement Mouhot for 
the initial discussion to look into the problem. 


( 3 ) 


2. Set up 

The stochastic differential equation 0 has the explicit solution 

X t =X 0 + i(l-e- At )K 0 + f t(1 — e~ x ^~ s ^)dW s . 
* Jo A 

V t = e~ xt V 0 + [ e- x{t ~ s) dW s , 

Jo 
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where Wt is the common Brownian motion. In this we separate the stochastic 
driving as ( A t ,B t ) given by the stochastic integrals 

A,-f 


(1 — e- A (*- s) )d W„ 


B t = f e- A(t - s) d W a , 
Jo 


which evolve over K with the common Brownian motion Wt- By Ito’s isometry 
(. A t ,B t ) is a Gaussian random variable with covariance matrix S (t) given by 


(4) 

(5) 

( 6 ) 


^ AA(t ) — 

S AB(t) = 


t~j( l-e- At ) + ^(l-e- 2At ) 


(l-e- At )--(l-e- 2A ‘) 


£bb(*) = - e 


-2A t ^ 


From this we calculate that the conditional distribution of A t given B t is a Gaussian 
with variance Yjaa{J) — 2 Yi 2 AB {t)Yi~j 3 1 B {t) and mean given by 

b) = EAs(0 s 5s(*) b - 


We write gA\B for the conditional density of A given B and gs for the marginal 
density of B. Hence 


(7) g{t,a,b) = g A \ B {t,a,b)g B (t,b) 

is the joint density of A and B. 

The last part of the set up is the change of variables we will need for the Mar¬ 
kovian coupling. We define new coordinates (Y, V) by taking the drift away 

f = x + If, 

V = V. 

The motivation for this change is the explicit formulas found in ([3| from which 
we see that Y is the limit as t -A oo of X t without additional noise. In the new 
variables, m becomes 

d Y t = id W t , 
d Vt = —\V t dt + d W t , 

for the common Brownian motion Wt- Note that the motion of Y t does not depend 
explicitly upon Vt and is a Brownian motion on the torus. 

It remains to show that these new coordinates define an equivalent norm on 
Txl. This follows from the triangle inequality and we have 

i) 

and the other direction is similar. Thus, the two norms are equivalent up to a 
constant factor that depends only on A. 


IX 1 - X 2 | t + IV 1 - H 2 | < IF 1 - F 2 | t + (l + 



3. Non-Markovian Coupling 

We wish to estimate how much the spatial variable will spread out over time. 
We will then use this to construct a coupling at a fixed time t which exploits the 
fact that a proportion of the spatial density is distributed uniformly. In order to 
do this we give a lemma on the spreading of a Gaussian density wrapped on the 
torus. 
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Lemma 7. For o > 0 consider the Gaussian density h on R given by 


h(x) = 


-\/27r<7 2 


0 -x*/2* 2 


and wrap it onto the torus T, i.e. define the density Qh on T by 

(9) (Qh)(x) = E h(x + 27 iLn). 


n£Z 


We have the following estimate on the spatial spreading 

P 


Qh(x) > 


where 


1-/3 = 


2t tL 

3e - CT 2 /2L 2 
2 g— < t 2 /2 L 2 ' 


Proof. By the definition of Q , the Fourier transform of Qh is for k € N given by 
(. FQh)(k) = f V h{x + 2nLn)e ikx / L dx 

n&L 

= f h{x)e ikx/L dx 

Jr 


= exp 


2 L 2 


where we have used the well-known Fourier transformation of a Gaussian. 
By the Fourier series we find that, for any ieT, we have 

1-/3 


Qh{x) - 


2ttL 2nL 


1 e -k 2 a 2 /2L 2 -ikx/L 


E 

ifci>i 


2t tL 


We want this to be positive. Therefore it is sufficient to show that 

a —k 2 a 2 / 2L 2 —ikx / L 


E* 

W>1 

We estimate the left hand side by 


E*- 

ifci>i 


<1-/3. 


< 2 e ~ k <y 2 /2L 2 = l_p 


k > 1 


where the final equality follows from summing the geometric series. 


□ 


We can now use this to construct a coupling at time t. We will use this coupling 
to prove exponential decrease in the Wasserstein distance. 

Lemma 8. Let fio, vq be probability distributions on TxR and let ((Xq, Vq ), (Xq, Vp)) 
be a coupling between them. Let t > 0 and /3 > 0 be such that for all b gR, 

( Q 9 A \ B ( t ,; b ))( a ) > 

where gA\B o,nd Q are defined by 0 and © respectively. Furthermore, let fit 
respectively be the distribution of the solution to the Fokker-Plank equation 0 
with initial data fi$ and vq respectively after time t. Then there exists a coupling 
((Xf, Vp), (Xf, Vfi)) between fit and v t satisfying 

E [(Vt 1 - V t 2 ) 2 } = e~ 2Xt E [(V 1 - Vo 2 ) 2 ] 















6 


HELGE DIETERT, JOSEPHINE EVANS, AND THOMAS HOLDING 


and 


E 


\Xl - X 2 | 2 ] < 2(1 - P)E 


l*o - X 2 | 2 ■ 


A 2 


(Vo 



Proof. Let us construct such a coupling. Split the distribution QgA\B as 
Q 9 A\B(t , a, b) = ^- + (1 - @)s(t, a, b). 

Then by assumption s is again a probability density for the variable a on the 
torus T. Let B be an independent random variable with density gB(t, b), let Z be 
an independent uniform random variable over [0,1] and let U be an independent 
uniform random variable over the torus. Finally let S' be a random variable with 
density only depending on B. 

With this define the random parts A 1 , A 2 of X 1 , X 2 as 


^=1 


A 2 =1 


z<[3 

Z</3 


u-xl 


u-x 2 


_L( l-e-^Vo 1 
l(l-e~ xt )V 2 


+ lz>pS, 


+ lz>/?S. 


By construction ( A 1 ,B) and ( A 2 ,B) both have law with density g(t,a,b) so that 
(X/, Vf 1 ) defined by 

Xl=Xl + ^{l-e- Xt )V£+A\ 

V 1 = e~ xt V£ + B, 


has law g t . and (X 2 , V 2 ) defined by 

X 2 =X 2 + j(l-e~ xt )V 2 + A 2 , 
V 2 = e~ xt V 0 2 + B, 


has law v t . 

Hence this is a valid coupling and we find 

E [(Vt 1 - V 2 ) 2 ] = e~ 2Xt E [(Ho 1 - H 0 2 ) 2 ] 

and 


E [I* 1 


X 2 | 2 ] = (1 - /3)E 



X 2 + -(1 - e-^ 1 


2 


X, 2 ) 


T 


and we can use Young’s inequality to find the claimed control. □ 

We now put these two lemmas together to prove Theorem [l], which states expo¬ 
nential convergence in the MKW W -2 distance. 

Proof of Theorem Q. Given any initial coupling of ((Xq, V^ 1 ), (Xq, H 0 2 )), we can use 
Lemma H to obtain a coupling ((X 1 , V^ 1 ), (X 2 , H t 2 )) of and v t . From explicitly 
calculating the variance of the distribution of A\B using Q, ©, ©, we see that 
the variance grows asymptotically as t/X 2 . Hence by Lemma [7] we can choose /? 
so that 1 — j3 0 exponentially fast with rate 1/2A 2 L 2 . This, combined with the 
control from the second lemma, shows that 

E [(Vt 1 - H t 2 ) 2 + IX 1 - X 2 | 2 ] < (V 2At + ce-V 2 ^ 2 ) E [(V 1 - V 2 ) 2 + \X% - X 2 | 2 ] . 

Taking the infimum over all possible couplings at time 0 gives the desired result. □ 


The explicit solution also allows to prove that the evolution is not a contraction 
semigroup. 
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Proof of Theorem We will prove the theorem by contradiction. Suppose 7 > 0 
and let a ^ b be two distinct points on the torus. Consider the initial measures 

MO = d x =a^v =0 


and 

UO = ^x—b^v— 0- 

Then the distance is W^Mo, ^ 0 ) = |a — 6|t- 

At time t the spatial distribution of Mt and 17 , interpreted in R, is a Gaussian 
with variance £aa which by the explicit formula Equation (0) can be bounded as 


Saa(^) A CAt 2 


for a constant Ca and t < 1 . 

Hence for d > 0 and t < 1 the spatial spreading is controlled as 
M(T\ [a - d,a + d\) x R) < 00 exp 

<C,Le x p(-C 2 t) 

for positive constants Ci and C 2 , where we have used the standard tail bound for 
the Gaussian distribution (see e.g. m Lemma 12.9]). 

For any d > 0 small enough that a ± d and b ± d do not wrap around the torus, 
any coupling between Mt and 17 must transfer at least the mass 

1 — m«((T \ [a — d, a + d]) xK)- v t {{T \ [b — d, b + d)) x R) 

between [a — d 1 a + d] and [b — d, b + d]. 

Hence the Wasserstein distance is bounded by 

> (|a - b\f - 2d) 2 ^1 - 20^—exp . 

Taking d = |a — b\ft 3 ^ 2 for t sufficiently small, this shows that 

> \a — 6 |t(1 - 2t 3/2 ) 2 ^1 - j^^-Vtexp . 

However, for all small enough positive t, we have 

(1 - 2 t 3 ’ 2 ) 2 > e -7t / 2 


and 


contradicting the assumed contraction. For the second estimate we use exp(— c/t) < 
(1 + c/t ) _1 = t/(c + f). □ 


4. Co-adapted couplings 

4.1. Existence. For Theorem^ we construct a reflection/synchronisation coupling 
using the drift-corrected positions Y t ‘. As the positions are on the torus we can use 
a reflection coupling until Yf and Yf agree. Afterwards, we use a synchronisation 
coupling which keeps Yf = l ) 2 and reduces the velocity distance. 

For a formal definition let ((Xg, V0, (Xq,Vq)) be a coupling between /i and v 
obtaining the MKW distance (the existence of such a coupling is a standard result, 
see e.g. in Theorem 4.1.]). For a Brownian motion Wf let (X 3 , Vf) be the strong 








8 


HELGE DIETEHT, JOSEPHINE EVANS, AND THOMAS HOLDING 


solution to m and define (A/ 2 , V t 2 ) as the strong solution with the reflected driving 
Brownian motion 


W 2 = | i 

t \ W} - 2W\ 


with the stopping time T = inf{i > 0 
analysis we introduce the notation 




t < T 
t > T. 

Y 2 } with Y'l from (fS|). For the 


M t = Yt - Y 2 

Zt = V? - V 2 . 


Then by the construction the evolution is given by 

(10) d M t = ^l t < T d W t \ 

(11) d Z t = -XZ t dt + 2 ■ l t <T&Wl, 


where M t evolves on the torus T. 

As a first step we introduce a bound for T. 

Lemma 9. The stopping time T satisfies 

(12) P(r>t|M„) = if; s T T exp(-l|^-()si» 

b —n \ / 


^ (2k + 1)|M 0 | t ^ 


Proof. As Mt evolves on the torus, T is the first exit time of a Brownian motion 
starting at Mq from the interval (0,27 tL). See [TO] (7.14-7.15)], from which the 
claim follow after rescaling to incorporate the 2/A factor. □ 


Remark 10. The second expression in (ED is obtained by solving the heat equation 
on [0, 2ttL\ with Dirichlet boundary conditions and initial condition Sm 0 ■ 

Lemma 11. There exists a constant C such that for any t > 0 the following holds 

(13) P(T > t\M 0 ) < C|M 0 | t (1 + r 1/2 )e- t/{ - 2>?L2) . 

Proof. Using ED and the inequality sin(a;) < x for x > 0, we have 


P(T > t\M 0 ) < -e-^ 2 ^ Y \ M °\i 2k+1 e -4#tm*L>) 
v 1 J - 7T ^ 2L 2k +1 

k=0 

< _|| Mo | T e-‘/( 2A2i2 ) (l + J™ e-* u2t /( 2x2L ^di 

= —\M 0 \ T e L } ^1 + y 8t/(X 2 L 2 )) 
<C|Af 0 | T (l + i" 1 / 2 )e- t/(2A2i2) 


where on the second line we have bounded the sum by an integral. 


□ 


Using these simple estimates, we now study the convergence rate of the coupling. 

Lemma 12. There exists a constants C such that for any t > 0 we have the bound 

( C|M 0 | T e- 2At 2A < 1/(2A 2 L 2 ) 

E [\M t \ 2 + \Z t \ 2 \(Z 0: M 0 )\ < \Z 0 \ 2 e~ 2Xt + l C\M 0 \ T (1 + t)e~ 2Xt 2A = 1/(2A 2 L 2 ) 

[c|M 0 | T e- t /( 2A2z ' 2 ) 2A > 1/(2A 2 L 2 ). 
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Proof. Without loss of generality we may assume that Zq and Mq are deterministic 
in order to avoid writing the conditional expectation. 

Applying Ito’s lemma, we find from m that 

d\Z t \ 2 = —2X\Z t \ 2 dt + 4 ■ l t < T Z t dWl + 2 • l t < T di. 

After taking expectations we see that 

(14) — E\ZA 2 = -2\mz t \ 2 + 2P(i < T). 

at 

By explicitly solving (ED and using Lemma 0 we obtain 

E\Z t \ 2 = \Z 0 \ 2 e~ 2Xt + 2e~ 2Xt [ e 2As P(s < T) ds 

Jo 

< \Z 0 \ 2 e- 2Xt + C\Mo\ T e- 2Xt 


—.it 

Let us bound I t . As the integrand is locally integrable, we have for a constant C 

It<c(^ + da 

Here the s -1 / 2 term can be bounded by 1 for s > 1 and for s < 1 the additional 
contribution can be absorbed into the constant. To bound the remaining integral 
we consider three cases: 

• 2A < 1/(2A 2 L 2 ): The integral (and I t ) are uniformly bounded, I t < C. 

• 2A = 1/(2A 2 L 2 ): The integrand is equal to 1 and I t < C( 1 + 1). 

• 2A > 1/(2A 2 L 2 ): The integrand grows and It < C{ 1 + e ( - 2A_1 /( 2A L ^*). 

In each case we multiply I t by e~ 2Xt to obtain the decay rate. In the first two cases 
this gives the dominant term with |Mo|t (as opposed to \Zq\) dependence, while 
in the last case it is lower order than the e~ t ^ 2X L ^ decay we obtain from E|M t || 
below. 

Next let us consider Using the finite diameter of the torus we have the 

simple estimate 

E\M t \l < tt 2 L 2 P(T > t). 

For t > 1 (say), we can use Lemma fl2l. to obtain 

E|M 4 |t < C|Mo| T e" t/(2A2i2) for t > 1 . 

This leaves when t < 1 where ED blows up. We instead use the martingale property 
of M t . Without loss of generality we may assume that Mq £ [ 0 , 7 tL]. Then as M t 
is stopped at T we know that M t £ [0, 27tL] for all t > 0. Hence, for any t > 0, 

E\M t \j < E\M t \ 2 < 2-KLEM t = 2i tLM 0 = 27rL|Af 0 |T 

by the martingale property. Combining the t < 1 and t > 1 estimates we have 

E\M t \ 2 < C'|M 0 | T e- 4/ ( 2A2z ' 2) for t > 0. 

This together with the bound for E|Z 4 | 2 provides the claimed bounds of the lemma 
and completes its proof. □ 

By the equivalence of the norms from (A, V) and (Y. V). this is the required 
coupling for Theorem^. 


[' e (2X-l/(2X*L*))s {1 + s -l/2 )ds 
JO 
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4.2. Optimality. In order to show Theorem 0, we focus on the drift-corrected 
positions Yf and Y 2 which behave like time-rescaled Brownian motion on the torus. 
For their quadratic distance we prove the following decay bound. 

Proposition 13. Suppose there exist functions a : (0,7 tL\ k > R + and f3 : [0,oo) H>■ 
R + with /3 £ T 1 ([0,oo)), such that, for any z £ (0,7tL] there exist two standard 
Brownian motions IF* 1 and IF* 2 on the torus T = R/(27 tLX) with respect to a 
common filtration such that |Wq — W 2 \ = z, and for t £ R + it holds that 

n\w} -w?$]<(a(z))*m 


Then with a constant c only depending on L, the function a satisfies the bound 

a(z) > c||/3||^i ( / [ o j0o)) \/i- 
From this Theorem 0 follows easily. 


Proof of Theorem^. Fix z £ (0,7rL] and consider the initial distributions p = 
<5x=o<V=o and v = Sx-zSv= o- Between p and v, there is only one coupling and 
W 2 (p,t/) = 2. 

If there exists a co-adapted coupling ((A* 1 , V) 1 ), (X* 2 , V 2 )) satisfying the bound, 
then Y ^/ x2 and Y 2 / x2 are Brownian motions on the torus with a common filtration. 
Moreover, 


r2|2 , + | F/-F * 2 ' 21 


n\Yt - Y?\l\<Cn\Xl - X?\ T - 

for a constant C only depending on A. Hence we can apply Proposition [13 to find 
the claimed lower bound for a. □ 


For the proof of Proposition [lil we first prove the following lemma. 

Lemma 14. Given two Brownian motions IF * 1 and IF * 1 on the torus with a common 
filtration, then there exists a numerical constant c such that 

n\Wl - IF* 2 | 2 ] > ce-^EOWo 1 - Wq It]- 

Proof. The natural (squared) metric \x — y\^ on the torus is not a global smooth 
function of x, y £ 1 as it takes x , y mod 2irL. Therefore we introduce the equivalent 
metric 

dj(x,y) = L 2 sin 2 , 

which is a smooth function of x,y £ R. Moreover, the constants of equivalence are 
independent of L , i.e. there exist numerical constants c\ and c 2 such that 

ci|a; - 2 /|t < d%(x, y) < c 2 \x - y\j. 


Now consider H* defined by 
H t = L sin 


IF* 1 - IF* 2 
2 L 


exp 


[IF 1 - IF 2 ]* 
4L 2 


As IF* 1 and IF 2 are Brownian motions, their quadratic variation is controlled as 
[IF 1 — IF 2 ]* < 4t. By Ito’s lemma 


dlL = - cos 


IF* 1 - IF 2 
2 L 


exp ^ W \ iL r ]t ) d(WF - IF 2 )*. 


Hence Ito’s isometry shows that 

/.* i 

E|H*| 2 =E|H 0 | 2 +E / -cos 2 

Jo 4 


IF* 1 - IF 2 
2 L 


exp 


[IF 1 - IF 2 ]* 
2 L 2 


d[W x - IF 2 ]* 


< E|H 0 | 2 +E J ex P {j^j dt 


< oo. 
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Therefore, H t is a true martingale and by Jensen’s inequality 

E[|-Ht| 2 ] > E[|£f 0 | 2 ]. 

Using the equivalence of two metrics, we thus find the required bound 

[w^-wV 


EHW * 1 - W 2 1 2 ] > c^E 


\H t \ exp - 


2 L 2 


> C 2 1 E [|iJ 0 | 2 ] exp 

> C 1 C 2 1 E[|W 0 1 - Wq |t] exp ■ 


□ 


With this we approach the final proof. 


Proof of Proposition Q3. Fix a £ (0,1), let z £ (0, nL\ be given, and by symmetry 
assume without loss of generality that Wq — W§ = IWq 1 — Wq\ = z. Then define 
the stopping time 

T = inf {t > 0 : W] - W 2 £ {az, nL)}. 


The distance can be directly bounded as 

E[|W/-W 2 |£]>P[T>t](a*) 2 . 


As /3 is integrable, it must decay along a subsequence of times and thus T must be 
almost surely finite. 

As Wf and Wf, considered on R, are continuous martingales, their difference is 
also a continuous martingale. By the construction of the stopping time, the stopped 
martingale (W 1 — W 2 ) t AT is bounded by nL and the optional stopping theorem 
implies 

P [W± - W| = nL] = Z ~ aZ . 


Since Brownian motions satisfy the strong Markov property, we find together 
with Lemma 14 


E 


/ OO 

|W/-W 2 | 2 dt 

POO 

> P [W± - W| = 7 tL]c {nL) 2 / e - 2t/i2 dt 

Jo 


z — az , ,9 Ij 

> — T — °W 2 — 

nL — az Z 

> CnZ 


for a constant C a only depending on a and L. 

On the other hand, integrating the assumed bound gives 

POO POO 

E / |Wl - W 2 \ 2 dt < (a(z)) 2 / 0(t)dt < (a(z)) 2 || / 9|Ui( [0 ,oc)). 

Jo Jo 

Hence 

C a z < (a(2;)) 2 ||/?|| L i( [ o, oo) ) 

which is the claimed result. □ 
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